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ENERGY  PARTITION  IN  THE  EXPLODING  WIRE  PHENOMENON 

ABSTRACT 

Streak  camera  and  oscillographic  circuit  damping  data  are  presented 
for  copper  vires  varying  In  diameter  from  5  to  8  ails.  A  maximum  of  specif io 
shock  wave  energy  in  the  induced  flow  is  found  at  a  wire  diameter  different 
from  that  of  a  minimum  in  the  total  damping  time  of  the  circuit.  This 
displacement  is  shown  to  be  caused  by  the  presence  of  residual  clroult 
resistance.  The  argument  is  based  on  a  critical  analysis  of  optimum,  damping 
conditions  in  the  exploding  wire  circuit.  A  maximum  of  apparent  energy  within 
the  contact  surface  appears  at  about  the  same  wire  diameter  as  the  minimum 
of  total  damping  time.  D1 5CUB  S  leu  Of  wll€  Implications  of  the  Taylor-  Lin 
similarity  theory  indicates  that  lack  of  similarity  of  the  flow  is  probably 
connected  with  the  displacement  of  the  maximum  energies  associated  with 
shock  wave  and  contact  surface. 
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15  INTRODUCTION 

We  continue  here  the  study  of  cylindrical  shock  waves  initiated  in  an 
1}* 

earlier  paper  .  There  it  was  shown  that  modification  of  the  rotating 
mirror  technique  makes  the  cylindrical  shock  wave  visible  over  a  considerable 
portion  of  its  patho  From  a  method  cjf  data  plotting  based  on  the  Taylor- 
Lin  similarity  theory  for  cylindrical,  blasts  one  finds  that  a  definite 
shock-energy  can  be  associated  with  each  trajectory* 

2) 

E.  David  '  has  recently  given  a  detailed  calculation  of  physical  pro¬ 
cesses  in  the  wire  explosion,  based  partly  on  the  experiments  of  4§£Lker'' 
and  partly  on  a  hypothetical  model  constructed  on  the  firmly  established 
principles  of  thermodynamics  and  electromagnetic  theory* 

We  are  not  concerned  here  with  the  finer  details  given  by  David,  but 
rather  with  the  cruder  task  to  establish  some  facts  relating  to  the  partition 
of  the  stored  condenser  energy  between  l)  energy  dissipated  in  the  induced 
fluid  motions  and  2)  energy  consumed  in  apparent  electrical  resistance  of 
the  circuit.  We  ahall  show  that  a  convincing  correlation  is  possible  between 
these  two  distinct  methods  of  energy  wastage,  with  the  result  that  the 
validity  of  our  method  of  assigning  a  shock  energy  is  thereby  rendered  more 
probable  * 


Superscript  numbers  denote  references  listed  at  back  of  the  report. 
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2. 


EXPERIMENTAL 


Fig.  1  shows  a  sequence  of  streak  camera  pictures  for  2  cm  copper  wires 
of  3,  4,  6.3  and  8  mils  diameter.  The  traces  were  obtained  by  use  of  the 
backlighting  technique  already  described  in  I. 

Condenser  energy  at  28  kv  is  118  Joules  and  the  ringing  frequency  of 
the  circuit  is  about  0.5  From  the  given  capacitance,  C  *  0.3  pf,  and 
the  dampi ng  curve  for  the  shorted  circuit  alone  (shown  in  I),  it  follows 
that  the  inductance  of  the  circuit  L  ip  about  l/3  ph  and  the  residual 
damping  resistance  R  is  approximately  l/4  ohm.  Damping  curves  for  each  wire 
are  shown  opposite  the  corresponding  flash  photograph  in  Fig.  1. 

The  shock  trajectory  is  most  clearly  seen  in  Fig.  1  (a)  where  separation 
from  the  contact  surface  occurs  by  about  0.8  p  sec.  For  the  4  mil  wire  of 

1  (c),  the  breadth  of  contact  surface  is  increased  and  separation  does  not 
occur  until  t  *  1.1  p  sec.  The  5  mil  wire  of  I  lies  intermediate  between 
this  case  and  the  6.3  mil  wire  of  1  (e),  where  separation  occurs  at  almost 

2  p  sec.  Here  the  shock  trajectory  is  less  distinct  then  in  the  previous 
two  cases  and  the  contact  surface  appears  to  have  gained  strength  relative 
to  the  shock.  Finally,  for  the  8  mil  wire  the  whole  phenomenon  is  narrower, 
appears  less  energetic  and  shows  a  visible  separation  of  the  shock  by  about 
t  *  0.6  p  sec. 

The  corresponding  oscillograms  show  interesting  detailed  behaviour, 
but  from  a  gross  point  of  view  demonstrate  mainly  that  total  damping  time 
passes  through  a  m-1  nl mirm  for  wire  diameter  close  to  6.3  mils.  Wires  of  9 
mil  diameter  and  higher  failed  to  give  a  luminous  effect. 

Because  of  the  interesting  behaviour  of  the  6.3  mil  wire,  the  complete 
flash  j.s  shown  in  Fig. 2.  Here  the  wedge  associated  with  the  interior 
facing,  second  shock  wave  discussed  in  I  is  very  prominent  and  occurs  at  t  = 

8  p  sec  well  after  the  electrical  oscillation  is  damped.  Just  ahead  of  the 
luminous  tip  of  the  wedge,  faint  traces  of  the  incoming  shock  may  be  seen. 
These  define  the  forward  nappe  of  a  cone  whose  outlines  may  be  detected  for 
perhaps  2  p  sec  prior  to  the  reflection  from  the  axis  represented  by  the 
wedge  tip. 
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3.  ANALYSIS  OF  DATA 


3.1  The  Energy  Plot 

Measurements  have  been  made  of  the  shock  and  contact  surface  trajectories 

shown  in  Fig.  1.  The  ordinate  in  these  flash  pictures  represents  radial 

distance  r  while  the  abscissa,  along  the  axis,  represents  elapsed  time  t. 

According  to  the  method  developed  in  I,  a  parabola- test  plot,  (2R)  vs 

(t-tn),  of  the  shock  data  should  yield  a  straight  line  if  the  experimental 

data  conforms  in  an  overall  fashion  to  the  Taylor-Lin  theory  for  an  ideal, 

k) 

cylindrical  blast  wave .  It  will  be  recalled  that  the  Lin  calculation 
is  based  on  the  assumptions  of  instantaneous  energy  release  along  a  mathe¬ 
matical  line  into  a  perfect  gas  of  constant  specific  heat  ratio  y  ■  1,4. 

The  actual  wire  explosion  releases  energy  into  a  finite,  massive  cylinder 
and  thence  into  an  imperfect  gas  (viz.  air)  whose  specific  heat  can  vary; 
nevertheless,  the  physical  situation  of  the  experiment  approximates  to  a 
reasonable  degree  the  assumptions  of  the  Taylor-Lin  theory. 

As  the  data  of  I  demonstrate,  the  cylindrical  shock  wave  formed  by  an 
exploding  wire  behaves  over  a  large  part  of  its  trajectory  like  an  ideal 
cylindrical  blast  wave  with  a  virtual  time  of  origin  tQ^o  and  an  apparent, 
axled,  energy  release  of  E  Joules  per  cm  easily  obtainable  from  the  slope 
of  the  data  plot » 

Rather  than  present  parabola- test  curves  for  each  size  of  wire,  a 

summary  plot  of  the  data  is  given  in  Fig.  3*  Here  the  quantity  Y  =  (2R)  /(t-to) 

has  been  plotted  for  both  the  shock  and  the  contact  surface  of  each  wire. 

As  may  be  seen  from  the  equation  for  the  shock  trajectory,  viz.,  1 

E  »  S  (j)  (E/p^1^  (t  -  t  J1/2,  a  plot  of  Y  vs  t  should  yield  a  horizontal 

straight  line  if  the  data  closely  approximates  an  ideal  blast  wave,  and 

2  l/2 

the  value  of  Y  will  be  given  by  m  =  4S  (E/pQ)  '  which  is  the  slope  one 
would  obtain  from  the  parabola-test  curve.  Assuming  values  of  S  [7)  and 
P0  appropriate  for  air,  we  can  readily  obtain  values  of  E  from  the  corresponding 
Y  values .  A  scale  of  energy  E  appears  at  the  right  of  Fig.  3. 
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Each  wire  size  is  represented  by  an  especial  symbol  designated  in  the 
legend  of  the  graph,  and  this  symbol  identifies  two  curves.  In  every  case 
the  higher  curve  represents  the  shock,  while  the  lower  represents  the  contact 
surface.  Horizontal  lines  have  been  drawn  to  fit  the  data  lying  at  or  beyond 
t  i'l  n  sec.  On  the  left  extremity  of  each  line  a  vertical  mark  indicates 
the  value  of  t  appropriate  to  the  fitted  data. 

The  points  early  in  time  tend  to  drop  from  higher  values  toward  the 
final  line.  This  transient  behavior  occurs  because  l)  the  early  data  of  a 
parabola- test  plot  define  a  knee- curve  whose  points  approach  the  approximating, 
inclined  straight  line  from  above,  and  2)  the  intercept  of  this  line  gives 
tQ7o ;  so  for  t  *  tQ  the  corresponding  Y  value  lies  at  positive  infinity. 

Values  of  Y  for  t  ?tQ  rapidly  approach  the  final  horizontal.  Because  this 
horizontal  line  defines  an  energy  value  for  the  blast  wave,  it  is  appropriate 
to  call  a  Y  vs  t  curve  an  "energy  plot". 

Several  interesting  facts  may  be  seen  by  inspection  of  Pig.  3.  In  the 
first  place,  Y  values  Indicating  highest  shock  wave  energies  of  approximately 
30+0.‘5  joules  /cm  occur  for  the  4  and  5  mil  wires;  but,  in  contrast,  the 
corresponding  contact  surface  curves  fall  well  below  that  of  the  6.3  mil  wire. 

The  4  and  5  mil  contact  surface  points  fall  steadily  further  below  the 
nominal  horizontal  line  drawn  at  14  Joules  /cm  as  time  increases  from 
t  =  1.5  (i  8ec.  The  3  mil  contact  surface  data  exhibits  the  same  tendency 
but  to  a  somewhat  less  pronounced  degree  beginning  to  deviate  toward  lower 
values  at  t  =  2  |i  sec. 

In  this  connection  notice  that  the  initial  times  for  the  contact  sur¬ 
faces  occur  earlier  than  those  for  the  shock  waves,  except  for  the  8  mil 
wire  where  the  two  coincide;  thus  the  virtual  starting  times  satisfy  the 
inequality  t  (contact  surface )^tQ  (shock  wave).  The  6.3  and  8  mil  wires 
are  notable  as  cases  where  both  shock  and  contact  surface  points  fall  well 
on  horizontal  lines  over  the  entire  range  of  data  obtained. 

3.2  The  Damping  Curves 

Inspection  of  the  oscillation  damping  curves  in  Fig.  1  and  that  of  for 
the  typical  5  mil  wire  given  in  I,  shows  a  further  distinction  to  be  made 
between  the  3,  4  and  5  mil  wires  when  compared  with  the  heavier  6.3  and  8  mil 
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wires .  For  the  first  three,  the  damping  curve  is  smooth  and  quite  continuous) 
for  the  latter  two  an  abrupt  interruption  occurs  after  which  the  oscillations 
are  almost  negligible. 


In  Table  1  we  collect  for  convenience  the  wire  resistance  at  20 UC . ,  the 
virtual  times  of  origin  of  contact  surface  (cs)  and  shock  (s),  the  total 
damping  uiiu8  f  and  the  energy  E  within  the  shock  wave. 

B 


TABLE  1 


fire  dlam. 
(mils ) 

■>  R_. 

(ohms) 

t  f  cs  ^  t  (ft} 

'o  ,  ~ 0  '  ' 

(n  sec) 

ncr 

time 

(4t) 

E 

1ofb.es /cm) 

3 

.075 

.06 

.51 

8.2 

19 

4 

.042 

.28 

.43 

7.2 

31 

5 

.027 

.33 

.46 

6.3 

30 

6.3 

.017 

.40 

.52 

3.5 

27. 

8 

.011 

.35 

.35 

3.5 

7 

9 

.0084 

4.0 

10 

.0068 

7.0 

The  damping  times  for  9  and  10  mil  wire  are  included  even  though  the  9  mil 
wire  produced  no  visible  flash  and  the  10  mil  wire  was  not  destroyed  by 
the  condenser  discharge. 


Inspection  of  the  lss"t  "two  columns  of  Table  1  elicits  two  facts: 
l)  a  maximum  in  the  indicated  energy  E-  occurs  for  a  wire  diameter  near 

B 

4  mils,  and  2)  a  minimum  of  toted  damping  time  occurs  for  a  diameter  between 
6.3  and  8  mils.  Furthermore,  the  damping=time  trough  is  50$  broader  than 
the  shock  energy  peak  on  a  scale  of  wire  diameter.  On  a  scale  of  diameter 
squared,  which  is  proportional  to  wire  mass,  the  trough  is  about  twice  the 
width  of  the  pfeak. 

3»3  Conclusions  from  the  Data 

The  data  just  discussed  indicate  certain  facts  relating  to  the  pro¬ 
duction  of  cylindrical  shock  waves  by  fine  copper  wires.  As  the  type  of 
wire  material  seems  to  play  no  q^sential  role,  providing  it  is  a  sufficiently 
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good  conductor,  these  facts  may  he  generalized  to  apply  to  Other  than  copper 
wires.  The  result  is  a  set  of  working  hypotheses  that  extend  our  view  of  the 
exploding  wire  phenomenon  and  provide  useful  tools  in  the  planning  of  future 
experiments.  Elevated  to  this  position,  they  may  be  stated  as  follows: 

Under  fixed  conditions  of  condenser  energy,  wire  material,  wire  length 
and  type  of  ambient  gas,  fluid  motions  produced  by  an  exploding  wire  have 
the  following  properties: 

1)  Both  shock  wave  and  contact  surface  follow  closely  parabolic  paths 
over  considerable  intervals  of  time. 

2)  A  maximum  of  Es,  the  specific  shock  wave  energy,  a  maximum  of  the 
Y  values  associated  with  the  contact  surface,  end  a  ml nimum  in 
damping  time  of  the  exploding  wire  circuit  occur  with  variation  of 
wire  diameter. 

3)  In  general  the  stationary  points  of  2)  do  not  necessarily  coincide, 
although  the  maximum  of  the  contact  surface  Y  function  appears 
close  to  the  minimum  of  total  damping  time. 


4.  DISCUSSION 

4.1  Maximum  Shock , Energy 

A  maximum  In  shock  energy  E-  considered  as  a  function  of  wire  diameter 
Is,  as  already  noted,  revealed  In  the  data  of  Table  1.  Its  presence,  although 
not  surprising  to  Intuition,  Is  not  easily  explained}  since  the  experimental 
concept  of  shock  strength  rests  here  on  interpreting  the  latus  rectum  of 
the  parabolic  streak  picture  in  terms  of  the  constants,  8(7)  (E/p^)1' \  given 
by  the  Lin  calculation.  The  goodness  of  approximation  of  the  actual  physical 
situation  to  the  hypotheses  of  the  ideal  fluid  calculation  given  by  Lin 
cannot  be  stated  with  precision  at  the  present. 

Accepting  the  interpretation  of  the  data  already  given  and  reading  from 
Table  1  the  maximum  energy  of  31  Joules/cm  it  is  apparent  that  the  2  cm  wire 
used  in  these  experiments  communicates  more  than  half  (2:;x  3l/ll8)  of  the 
stored  condenser  energy  into  the  fluid  motion  associated  with  the  shock  wave. 
End  defects  have  been  purposely  neglected  in  arriving  at  this  eatimatei 

4.2  Circuit  Resistance  and  Minimum  Damping  Time 

The  presence  of  a  minimum  damping  time  at  about  7  mils  wire  diameter 

is  readily  seen  from  Table  1.  This  fact  implies  that  the  resistance  of  the 

7  mil  wire,  R  ,  plus  the  residual  resistance  of  the  circuit  R  ,  forms  a 
w  c 

total  that  is  in  some  sense,  as  yet  undefined,  well  matched  to  the  L-C 
parameters  of  the  circuit. 

As  already  mentioned  in  I,  the  residual  resistance  Rq  is  about  l/4  fl. 

Values  of  R  with  the  exploding  wire  shorted  out,  calculated  from  the 
c 

logarithmic  decrement  of  successive  current  peaks,  show  R  initially  to 

c 

start  at  about  0.25  fl  and  to  increase  slightly  up  to  0.35  fi  at  the  tail  of. 
the  damping  curve. 

A  similar  calculation  made  from  the  damping  ^curves  of  the  3,  4  and  5 
mil  wires  shows  Rc  starting  at  about  0.30  0  and  increasing  to  at  most  0.45  fl; 
which  evidence  may  be  interpreted  to  mean  that  by  the  time  a  clear  oscillo¬ 
graphic  record  is  available  for  these  wire  sizes,  that  is  after  about  1  p 
sec  has  elapsed,  an  electric  discharge  has  been  established  in  the  expanding 
metal  vapor  and  the  subsequent  damping  is  dominated  by  this  discharge  and 
the  resistance  residing  in  the  remainder  of  the  circuit. 


11 


2) 

David  estimates  that  the  conductivity  of  the  copper  wire  should  drop 
suddenly  after  the  diameter  has  increased  by  a  factor  of  two  to  four.  When 
this  condition  is  reached  metallic  conduction  should  no  longer  exist,  but 
electron  transfer  could  occur  perhaps  only  by  tunnel  effect  modified  by 
local  fields.  After  the  expansion  has  proceeded  sufficiently  far,  the 
pressure  of  the  metal  vapor  will  drop  to  a  point  where  arc  discharge  can 
occur  if  enough  energy  remains  in  the  circuit  and  the  conductivity  will 
increase  again  to  a  high  value.  In  the  event  that  little  energy  remains 

9 

by  virtue  of  its  rapid  expenditure  in  i“R  heating,  no  arc  can  form  after 
the  wire  diameter  exceeds  2-1*.  times  its  original  value.  The  effect  will  be 
like  that  of  opening  a  switch. 

The  oscillograms  for  the  6.5  and  8  mil  wire  indicate  -Just  such  an  effect 
after  about  1  and  2  p,  sec  respectively.  The  damping  records  show  that  these 
wires  expand  the  condenser  energy  more  rapidly  than  the  5»-  4,  5  ail  series. 

We  infer  the  cause  to  be  better  impedance  match  of  the  exploding  wire  to 
the  circuit.  The  streak  picture  for  6.5  ails  given  in  Pig.  2  tends  to  bear 
out  this  interpretation  if  the  contact  surface  and  second  shock  are  regarded 
as  being  more  energetic  than  in  the  other  cases. 

4.5  Choice  of  Optimum  Circuit  Resistance 

The  evidence  just  discussed  supports  the  assumption  that  the  6.5  *  9 
mil  wires  provide  nearly  optimum  electrical  resistance  for  maximum  power 
and  energy  consumption  in  the  circuit.  The  effective  resistance  is  clearly 
far  below  the  critical  damping  resistance j  for,  even  though  damped  to  a 
minimum  total  time,  the  circuit  is  still  in  the  oscillatory  regime. 


In  order  to  carry  the  discussion  further,  criteria  bearing  on  the 
question  of  maximum  power  dissipation  in  an  R-L-C  circuit  are  needed. 
Investigation  of  the  available  literature  shows  that  even  for  the  case  of 
fixed  circuit  resistance  R,  the  conditions  for  maximum  power  dissipation 
have  not  hitherto  been  precisely  defined. 


-  — 


xn  connecT/iun  vi  ui  a  study  of  exploding  wires  made  in  1926,  Anderson 


and  Smith'*) 

circuit  with  constant  parameters  to  be  approximately  R  =  (L/C) 


give  the  condition  for  maximum  power  dissipation  in  an  R-L-C 

1/2 
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Because  they  provide  no  discussion  of  this  result  and  because  their  subsequent 
statements  about  maximum  current  are  open  to  objection,  the  question  has  been 
re-examined  and  some  additional  results  obtained.  The  more  precise  and 
extended  analysis  is  given  below  in  Appendix  I. 

There  it  is  found  that  for  maximum  power  dissipation  at  the  first  current 
peak  the  total  circuit  resistance  should  bdj  R  =  1.10  (L/C)1/2.  If  the  teat 
losses  are  i  enough  so  that  the  process  can  be  considered  adiabatic, 
then  energy  dissipation  in  a  specified  interval  may  be  of  primary  importance. 
For  example,  if  this  interval  extends  past  the  first  current  peak  to  terminate 
when  the  sine  factor  in  the  current  expression  is  unity,  then  R  =  1.24  (l/C)1^ 
should  i.bold.  If  a  longer  interval  is  important,  the  appropriate  numerical 
factor  will  be  larger  than  1.24.  As  the  interval  becomes  infinite,  R 

l/p 

approaches  the  critical  damping  resistance  given  by  2  (L/C)  '  . 

Now  it  is  obvious  that  the  resistance  of  an  exploding  wire  is  not  con¬ 
stant;  so  that  the  criteria  given  above  cannot  strictly  apply  except  in  the 
case  that  power  consumption  at  maximum  current  is  the  dominating  term.  Even 
here  no  energy  can  be  dissipated  save  during  a  finite  interval.  Therefore 
we  conclude  that  the  optimum  average  resistance  of  the  exploding  yite  circuit 
should  satisfy  the  inequality  1.10£Ropt  (l/C)”1/2^  1.30  where  the  lower  bound 
is  that  set  by  maximum  power  dissipation  at  the  first  current  peak  and  the 
upper  bound  is  fiet  slightly  to  exceed  the  value  for  maximum  energy  dissipation 
in  the  approximately  quarter-cycle  interval  already  discussed.  The  oscillogram 
for  the  6.3  mil  wire  suggests  an  interval  before  current  interruption  some¬ 
what  larger  than  a  quarter  cycle  innsupport  of  a  choice  of  upper  bound  greater 
than  1.24. 

With  the  L,  C  values  for  the  present  circuit  this  inequality  works  out 
to  be  l.l6&Ropt£l.36  ohms.  As  already  seen  the  residual  circuit  resistance 
is  approximately  l/4  -  l/3  fij  so  by  difference  the  exploding  wire  should 
furnish  0.9  -  1.2fiO  in  order  that  the  total  resistance  most  rapidly 
dissipate  the  stored  condenser  energy. 

2) 

Relying  here  on  David’s  discussion  '  of  the  resistance  changes  induced 

by  heating  of  the  copper  wire,  we  estimate  that  wire  resistance  may  increase 
2  “ 
by  a  factor  of  10  for  a  rise  from  room  temperature  up  to  4000°K.  Applying 
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this  factor  to  the  cold  resistance  given  in  Table  1  wa:  find  that  the  6.3,  7 
and  8  mil  wires  would  have  high  temperature  resistances  in  or  near  the  range 
of  values  suggested  as  optimum  in  the  previous  paragraph.  Wires  of  5  mil 
diameter  or  smaller  would  appear  to  have  high  temperature  resistances  too 
high  by  a  factor  of  two  or  more  and  larger  than  the  critical  damping  resistance. 

One  would  not  expect  the  trough  surrounding  the  minimum  of  damping  time 
to  be  narrow  since  the  maximum  of  power  dissipation  at  the  first  current 
peak  will  merge  smoothly  with  the  relative  maxima  of  energy  dissipation  as 
the  effective  time  interval  increases  from  zero.  (See  discussion  in  Section 
2  of  Appendix  I.)  The  experimental  damping  times  are  in  agreement  with 
this  qualitative  prediction  from  theory;  for  damping  times  from  6.3  to  9  mils 
are  roughly  the  same  and  the  minimum  is  broad  compared  with  the  maximum  of 
shock  energy. 

It  appears  highly  probable  that  the  nearly  optimum  wires  reach  temperatures 
of  several  thousand  degrees  Kelvin,  perhaps  even  the  critical  point  for 
copper  under  influence  of  magnetic  pinch  pressures,  and  that  their  high  temper¬ 
ature  resistance  values  are  close  to  those  demanded  for  maximum  power  or 
energy  dissipation.  The  data  appear  to  favor  values  higher  than  the  lower 
hound. 

Effect  of  Residual  Circuit  Resistance 

We  develop  here  an  argument  to  show  that  the  displacement  of  the  shock 
energy  peak  from  the  minimum  of  damping  time  is  entirely  caused  by  the 
residual  resistance  in  the  circuit.  Conceding  temporarily  the  truth  of  this 
proposition  then  if  this  resistance  is  made  to  vanish,  one  would  expect  the 
shock  with  greatest  axial  energy  release  and  the  most  rapid  dissipation  of 
condenser  energy  to  occur  at  the  same  wire  diameter,  a  criterion  which  will 
clearly  be  of  importance  where  luminous  effects  (e.g.  light  sources)  of 
minimum  duration  are  desired. 

As  may  be  seen  from  Appendix  I  the  necessary  and  sufficient  conditions 

that  the  A  power  in  the  circuit  be  stationary,  i.e.  dP  *»  0,  are  that 

Si/St  ■  0  and  Sin  i/SR  -  -l/2R.  As  before  the  resistance  of  the  circuit 

is  the  sum  R  *  R  +  R  where  R  denote  residual  and  exploding  wire  resistance 
c  w  c 
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respectively.  Wire  resistance  R  is  directly  proportional  to  length  1  and 
inversely  proportional  to  radius  r  squared;  thus  R  =  kl/r  .  The  power 

v  g  2 

expended  Per  unit  length  of  wire  is  given  by  pw  =  i**/1  *  k  (i/r)  . 

Intuitively  it  is  obvious  that  p^At  must  be  related  to  the  specific  shock 
energy  E  of  the  Lin  theory  because  the  energy  expended  in  i2Rw  heating 
partially  reappears  to  supply  the  energy  of  the  fluid  motions  induced  by 
the  explosion.  Here  At  represents  a  small  time  interval  about  the  first 
current  peak. 

By  logarithmic  differentiation  weifind  d  In  p^  ■  2  d  In  (i/r);  but  at 
the  instant  of  power  maximum,  current  is  a  function  of  R  alone  and  in  particu¬ 
lar  d  In  i  *  -(l/2)  d  In  R.  Combining  these  two  results  leads  us  to 

2 

d  In  p  *  -d  In  (Rr  ).  A  stationary  point  for  py  /  0  can  only  exist  when 
the  right  hand  side  of  this  equation  vanishes.  Otherwise,  a  ncn=zerc 
value  will  obtain  and  an  incremental  change  in  wire  diameter  (resistance) 
can  alwygrs  be  found  which  will  increase  ptt. 

2  2 

Prom  the  definitions,  Rr  =  RQr  +  kl  and  we  have  finally,  dp^/p^  => 

-  (2R^/R)  (dr/r),  from  which  it  is  obvious  that  dp__/dr  £0  and  vanishes  if 
Rq  *  0.  We  seen  then,  as  stated  at  the  beginning  of  this  section,  that 
non-vanishing  residual  circuit  resistance  implies  a  separation  between  the 
minimum  of  damping  time  and  EwRy  for  the  shock  wave;  furthermore,  decreasing 
the  wire  diameter  from  the  value  necessary  to  satisfy  dP  »  0  causes  pw  to 
increase . 


Assuming  from  the  data  of  Table  1  that  E  occurs  at  k  mils  and  mini- 
mum  damping  time  at  6.3  mils,  with  R  =  l/l*  (1  and  R  ®  1  fl,  we  find 
dp„/p„  “  2  x  (l/5)  x  (2. 3/6. 3)  which  works  out  to  a  15#  increase  in  p._ 
upon  decreasing  the  wire  diameter  from  6.3  to  k  mils.  Inspection  of  Fig.  3 
shows  a  corresponding  increase  in  the  shock  wave  energy  of  about  12#;  which 
in  view  of  the  uncertainties  in  the  values  used,  may  be  considered  as 
remarkable  agreement. 


It  seems  clear  then  that  the  6.3  mil  wire  corresponds  closely  to  a 
damping  situation  iln  which  maximum  energy  is  dissipated  during  an  interval 
about  the  current  peak.  The  fact  that  the  neighboring  wires  that  show 
efficient  short  damping  times  have  lower  cold  resistances  and  presumably 
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lower  high- temperature  resistances  suggests  that  the  maximum  Integrated 
power  dissipation  of  the  upper  bound  may  apply  rather  than  the  maxi  mum  power 
at  peak  current.  Nevertheless  the  foregoing  analysis  based  on  the  latter 
matching  criterion  brings  out  excellent  agreement  between  the  predicted 
specific  power  increase  in  the  wire  and  the  observed  energy  increase  in  the 
shock  wave. 

Because  this  agreement  is  found  between  quantities  obtained  by  the  two 
very  different  approaches  taken  through  l)  the  electrical  theory  of  the 
R-L-C  circuit  and  2)  fluid  mechanical  theory  of  the  cylindrical  blast  wave, 
our  confidence  is  strengthened  in  the  use  of  the  second  approach  to  derive 
energy  values  for  the  experimental  shock  waves. 

It-.  5  Contact  Surface 

We  revert  now  to  the  question  concerning  the  connection  between  an 
apparent  energy  within  the  contact  surface  and  lack  of  similarity  in  the  flow. 

The  similarity  theory  assumes  that  the  cylindrical  shock  wave  is  a 
constant  energy  surface  of  the  flow.  ^Examination  of. the  energy  integral 
[Lin's  Eq.  (23)]  shows  that  for  arbitrary  radius  r  interior  to  the  shock, 
the -energy  contained  within  r  is  a  function  only  of  t)  -  r/R.  Thus  parabolae 
given  by  tj  *  const. £  1  are  curves  of  constant  energy  in  the  r-t  plane. 

On  the  other  hand,  particle  paths  in  the  r-t  plane  define  curves  of 
constant  mass,  since  by  definition  these  paths  may  not  intersect.  Because 
each  particle  experiences  a  radial  deceleration  after  passage  of  the  shock, 
each  inner  particle  must  be  doing  work  on  the  next  outer  particle  of  fluid. 
Therefore  energy  1b  flowing  outward  faster  than  the  particles  on  any  given 
cylinder  of  radius  r,  and  particle  paths  cut  constant  energy  parabolae  in 
the  r-t  plane  in  the  direction  of  decreasing  tj  values  as  time  increases. 

Now  the  contact  surface  experimentally  determined  from  streak  camera 
pictures  may  reasonably  be  assumed  to  approximate  the  outer,  limiting  particle 
path  of  the  luminosity.  If  the  flow  is  truly  similar  in  the  Taylor-Lin 
sense,  then  the  contact  surface  should  not  coincide  with  a  constant  energy 
surface  given  by  ^  *  const.,  but  should  show  a  trend  of  tj  toward  decreasing 
values.  In  contradiction  to  this  conclusion  the  6.3  and  8  mil  data  of  Fig.  3 
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appear  to  lie  well  on  curves  tj  a  const..  (The  data  for  the  smaller  wires 
show  a  decreasing  trend  of  tj  with  time.)  Thus  we  may  certainly  say  that 
the  flows  generated  by  the  larger  wires  fail  to  exhibit  similarity. 

On  the  other  hand,  these  flows  have  contact  surfaces  which  plot 
like  one  of  the  constant  energy  surfaces  of  the  ideal,  similar  flow 
corresponding  to  the  experimental  shock  wave.  Thus  by  comparison  between 
the  experimental  flow  and  the  ideal,  similar  flow  of  equal  shock  energy, 
the  contact  surface  may  be  assigned  an  energy  value  which  is  the  same  as 
the  energy  inside  the  tj  *  const,  surface  of  the  ideal  flow  on  which  the 
experimental  contact  surface  points  fall.  While  this  energy  may  not  represent 
accurately  the  total  energy  residing  within  the  experimental  contact  surface, 
it  is  nevertheless  a  well  defined  quantity  derivable  from  the  contact 
surface  trajectory,  and  may  be  valuable  for  comparison  of  one  contact  surface 
with  another.  We  adopt  the  foregoing  definition  of  energy  within  the  contact 
surface  for  all  subsequent  discussion. 

To  examine  the  contact  surface  data  for  the  3,  k  and  5  mil  wires  more 
closely,  we  need  detailed  information  about  the  particle  trajectories  of 
the  similarity  flows.  To  this  end  recall  that  for  the  similar  solutions 
fbundu  by  Lin1'',  the  velocity  of  a  particle  at  radius  r  is  given  by  u  »  0(tj)  U, 
where  U  *  dR/dt  «=  R  gives  the  velocity  of  the  shock  and  0(  t|)  is  a  numerically 
determined  function. 

The  particle  paths  may  be  Obtained  by  integrating  the  Lagrangian 

•  » 

equation  of  motion  r  =  u  =  0R.  With  the  help  of  r  ■=  tj$  +  fjR  we  find  the 

first  order  differential  equation  dT)  »  (0  -  tj)  d  In  R.  From  Lin!s  Table  I 

we  t'f ind  that  (0  -  tj)  is  always  negative.  As  a  result  we  have  dr]/dt  * 
U(dT}/dR)^L0  and  the  decreasing  trend  of  tj  already  noted  is  verified  again. 

The  fractional  change  in  tj  is  found  from  drj/l  ®  tj  "‘"(0  -  tj)  (dR/R) 
from  which  it  follows  with  the  help  of  Lin's  tabulated  values  that 

.275  (dR/R).  The  fractional  fahange  in  tj  divided  by  that  of  R 
varies  from  0.16?  at  tj  *  1  to  O.275  at  tj  =  0.6  and  appears  to  remain 
constant  at  this  value  as  tj-*0.  Thus  we  see  that  tj  for  a  particle  trajectory 

is  a  slowly  varying  function  of  R  and  is  to  this  extent  crudely  approximated 

by  setting  tj  =  const.,  as  was  suggested  in  I. 
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The  particle  paths  can  be  obtained  by  a  numerical  integration  of  the 
equation  (0  -  x)'1**  -  in  (R/Rx)  -  (l/2)ln  [(t  -  tQ)/(t;L  -  tQ)J  , 

where  R7  and  t7  are  the  shock  radius  and  time  respectively  at  which  the 
shock  passes  over  the  particle  in  question,  and  tQ  is  the  virtual  time  of 
origin  of  the  shock  as  before.  Values  of  the  Integral  In  (r/r..  )  have 
been  obtained  by  Stepson's  rule  from  the  numerical  results  in  Lin's  paper4 
and  are  presented  in  Table  2.  Using  these  values  and  the  experimental 
values  for  t^  and  tQ,  we  have  calculated  theoretical  and  experimental  values 
for  the  particle  trajectory  defined  by  the  constant  surface  data  for  the 
and  5  mil  wires.  These  data  are  presented  in  the  form  qf  ratios  ] 
r](exp.)/q  (calc.)  versus  time  in  Table  3.  Of  the  three  wires  the  L  mil  comes 
nearest  to  similarity  being  only  2$  low  beyond  1  |x  sec.  The  5  mil  wire  is 
close  to  similarity  in  the  latter  part  of  the  range  while  the  3  nil  wire 
is  never  closer  than  8$  and  as  far  as  16$  below  the  calculated  values. 

The  experimental  trajectories  nearly  always  lie  inside  those  calculated 
from  theory.  Note  that  the  nearest  approach  to  a  similar  flow  occurs  for 
the  wire  with  the  most  energetic  shock  wave,  viz.,  the  ^  mil  wire. 


Turning  now  from  discussion  of  the  lack  of  similarity  in  the  flows, 


We  conclude  by  eXealxiicg  briefly  the  xu^B^dn  3f  the  relative  aSCUnt  Of 


energy  within  the  contact  surface.  The  definition  given  above;  which 
relates  the  experimental  contact  surface  to  an  T]  *  const,  curve  of  the  ideal, 
similarity  flow  of  equal  shock  energy,  provides  a  definite  number  for  the 


contact  surface  but  no  assurance  that  this  number  represents  the  real 
energy  within  that  surf  ewe. 


The  fact  that  the  highest  value  occurs  for  the  6.3  mil  wire  whose 
streak  picture,  Fig.  2,  shows  Vhat  appears  to  be  the  most  energetic  phenomena 
within  the  contact  surface,  encourages  the  hypothesis  that  there  is  a 
proportionality,  if  not  an  identity,  between  the  actual  energy  within  the 
contact  surface  and  the  graphically  determined  value  obtained  from  the 
energy  plot. 


Since  the  mar-t of  shock  wave  and  contact  surface  Y  values  do  not 
occur  at  the  same  wire  diameter  we  are  forced  to  suppose  that  the  apparent 
energy  within  the  contact  surface  is  regulated  not  only  by  phenomena 
responsible  for  the  shock  wave,  but  also  by  other  factors  at  present  im¬ 
perfectly  known,  which  have  to  do  with  the  electrical  and  mechanical  phenomena 
of  the  wire  explosion. 


18 


The  experimental  data  of  Figs.  1,  2  and  3  seem  to  support  these  hypotheses ; 
for  the  most  quickly  damped  condenser  discharge  iB  associated  with  the 
contact  surface  of  highest  indicated  energy  and  with  the  most  energetic 
appearing  second  shock  wave.  The  most  energetic  shock  wave,  on  the  other 
hand,  is  associated  with  a  contact  surface  definitely  below  the  contact 
surface  maximum  on  the  energy  plot.  This  suggests  that  while  the  maximum 
E  value  for  a  shock  wave  appears  at  a  smaller  wire  diameter  k)0CcfcUJ3G  0±  the 
distribution  of  electrical  energy  between  internal  and  wire  resistance, 
the  condition  of  minimum  damping  time  occurs  at  the  same  wire  size  as  flow 
phenomena  of  maximum  kinetic  energy  within  the  contact  surface.  That  is  to 
say,  maximum  dissipation  of  electrical  power  is  accomplished  partly  through 
the  simultaneous  appearance  of  fluid  dynamic  effects  different  from  those 
responsible  for  the  generation  of  the  shock  wave.  Not  much  more  can  be  said 
at  this  point;  for  clearly,  further  experimental  and  theoretical  information 
is  needed. 
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TABLE  2 


Table  2  Numberical  values  of  the  Integral  J?.-  (0  -  x)  ^dx  versus 

the  parameter  tj  *  r/R. 


n 

J2 

n 

X 

7]  X 

n 

X 

1.000 

0 

•  94 

.332 

.82  .924 

•55 

2.420 

•  995 

.029 

.93 

.384 

.80  1.021 

.50 

2.766 

-990 

.059 

.92 

.435 

.78  1.119 

.40 

3.569 

.985 

.088 

•  yj- 

•  485 

.76  1.218 

•  ^ 

4.627 

.980 

.116 

.90 

•  535 

•74  1.319 

•  20 

6.097 

•  97 

.172 

.88 

.633 

.70  1.539 

.10 

8.673 

.96 

.226 

.86 

.731 

.65  1.811 

•  95 

.280 

.84 

.827 

.60  2.103 

TABLE  3 

Comparison  of  calculated  and  experimental  trajectories  for 
the  contact  surfaces  of  the  3*  1  and  5  mil  wires. 

t  (u  sec) 

.003" 

TKexp;  /  rj^caxc.) 

.004"  .005” 

1.0 

.88 

.93 

•  93 

1.5 

.89 

.98 

•  93 

2.0 

.92 

•  99 

ni. 

•  7*+ 

2.5 

•  89 

.98 

•97 

3.0 

.86 

.98 

.98 

3.5 

.84 

- 

1.01 
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a.  .003"Copptr  Wir«  o. 


Fig.  1. 


Rotating-mirror  pictures  and  damping  curves  for  sequence  of  copper 
wires  under  constant  ambient  conditions  and  condenser  energy.  The 


M-  sec  scale  in  (g)  applies  to  all  flashes. 


jiacn  norizoniai  division 


of  the  oscilloscope  scale  corresponds  to  2  q  sec. 
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Fig.  2 


.  Rotating-mirror  picture  for  exploded  copper  wire  of  6.5  mil  diameter  to 
show  details  of  the  wedge  luminosity  representing  the  second  shock  wave. 


2k 


t  fJL  sec 

fig.  .5-  Energy  plot  of  shock  wave  and  contact  surface  data  for  the  sequence  of 
exploded  copper  wires  shown  in  Fig.  l. 


APPENDIX  I 

POWER  DISSIPATION  IN  AN  EXPLODING  WERE  CIRCUIT 


1.  MAXIMUM  AT  PEAK  CURRENT 


Anderson  and  Smith^  discuss  the  typical  R-L-C  circuit  used  to  produce 
the  heavy  current  pulse  under  which  a  thin  metal  wire  explodes.  They  state 
without  proof  that  a  maximum  in  the  power  dissipated  by  the  circuit  will 
occur  when  the  dieeuit  resistance  R  has  a  value  near  (l/C)^^.  It  may  be 
inferred  from  their  analysis  that  current  i  is  to  be  maximized  simultaneously 
with  power  although  in  the  later  step  of  obtaining  the  maximum  squared 


uaptlon  is  ignored c 


icao wuo , 


An  explicit  solution  to  this  same  problem  is  given  here  for  the  following 
L)  a  complete  discussion  is  not  given  by  Anderson  and  Smith  and  cannot 
readily  be  found  elsewhere  in  the  literature,  2)  the  equation  given  by 
Anderson  and  Smith  for  the  maximum  current  squared  is  not  precise  for  reasons 
pointed  out  belowj  thus  sdme  doubt  attaches  to  their  earlier  assertion  about 
the  best  value  of  R,  and  3)  consideration  of  integrated  power  leads  to  some 
new  "best  values"  of  R  which  may  sometimes  be  preferred. 

2  2  l/2 

Defining  the  quantities  u  =  R/2L  and  to  *  (l/LC  -  R  /tL  )  '  01 

write  for  the  current  in  the  circuit 


/,» /  \  -at  , 

e  sin  ait 


( i  > 
^ ) 


L  uc  d  Ou  UlviiCU  Jr  points  of  i,  or  In  i  for  convenience,  are  determined  from 
d(ln  i)/dt  =  to  cos  cot/sin  cot  -  a  =  0  ,  (2) 

which  gives 


cot  cot  =  a/co 
From  (3)  the  sequence  of  times  I 


(3) 


maxima  and  minima  of  current. 


tA  J-  may 


be  determined  to  locate  successive 


Instantaneous  power  dissipated  in  the  circuit  is  A  and  the  condition 
d(lT0/SR  »  0  implies  that 

d(ln  i)/dR  =  -  l/2R  .  (1+) 
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With  (l)  and  the  definitions,  Eq.  (4)  leads  to  the  relation 
(t  cot  out  -  i/cd ) ( Scd/SR )  -  t(da/&R)  +  l/2R  =0  , 

or.  with  the  nartials  explicitly  evaluated,  to 


(5) 


l/2a  -  (a/co)(t  cot  cot  -  l/co)  -  t  =  0  .  (6) 

Enns.  ( and  (6)  determine  the  conditions  necessary  to  maximize  power  with 

_  x - -  \  f  >  —  /  "  •»  ■  -  — 

respect  to  i  and  R  simultaneously,  Eqn.  (3)  is  itself*  necessary  to  maximize 
i  alone.  We  consider  only  the  first  maximum,  since  the  wire  is  gone  shortly 
thereafter . 

If  (3)  he  solved  for  t  and  used  in  (6)  with  definition  p  =  o>/a,  there 
results  the  transcendental  equation 

tan_1p  =  (p/2)(p2  +  2)(p2  +  l)”1  .  (7) 


1 _ j.  j _ 


-  _ _ 


axaminaoion  or  (7)  together  with  (3)  scows  xnax  mere  is  an  inrinity  ox 
increasing,  positive  roots,  beginning  with  zero,  which  corresponds  to  the 
series  of  successive  stationary  current  values.  A  graphical  solution  gives 

■n  •  1  Cl  C  +U  A  -P  4  MAM  KAMA  MAA+  PmAM  +V4  A  4  +  1  At*A  +W  A  + 

—  x«  yxy  xwju  wuc  x  -1.x  o  u  iiuu-cciu  x  ww  \j  o  i'iuiu  uixxo  xu  x  uxxun o  uuctu 

R  =  2(1  +  p2)-1/2  (L/C)1/2  =  1.10  (L/C)1/2.  If  one  defines  2(1  +  p2)-1/2  by 

f(p)  in  the  general  case,  then  0  £  f  2  comprises  all  values  for  which 
oscillations  may  occur  in  the  circuit.  Values  f  =  0,  2  correspond  to  the 

luiucuupcu  cuiu  un  biuanjf  aouijjcu.  iui  wuiuu  y  ‘•W;  \j  rco^cwoiYCiJ,  cum 

f  =  1.10  represents  the  resistance  in  which  power  is  dissipated  most  rapidly 
at  the  first  current  maximum.  It  is  now  clear  that  Anderson  and  Smith's 

rwA  rr4  m  o  1  occoy«Hnn  4  a  omH  wol  on+  4-  r\  covl  nrr  4'Vi  o  +  mov^nmni  t**  rYr.ro  v*  ^4  cc-It\d-H  r\y~, 

WX  X^AXiMO.  WM  WS.J.  WXWU  J.  U  WV^M4.  VV  UMJ  VilUV  A  VA  UUAAXiUVUU  ^/vn^j.  UXUOX^/UUXWi^ 

at  the  first  maximum  current,  f  must  be  close  to  one;  furthermore  the  explicit 
solution  shows  that  f  is  actually  10$  greater  than  one.  For  the  successive 
current  maxima  decreases  monotonl call v  to  zero. 


The  expression  for  current  may  be  rewritten  as 

i  =  (VapL)  e^^/p  sin  cot  .  (8) 

The  stationary  values  occur  when  (3)  is  satisfied,  i.e.  when  oyt_  =  tan_1p  +:  nut 

III 

with  m  =  0,  1,  2...  and  sin  aytm  =  (-l)m  p  (l  +  p2)-1/2.  With  these  relations, 
(8)  becomes 
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,  .  «m,_-  /  .  x  /.  2\-l/2 

im  =  (-1)  ^V/qcl;(i  +  p  ;  '  exp 


-p  itui  p  -  uoi 


/  n  \ 
\y ) 


which  furnishes  the  first  maximum  when  m  =  0.  From  the  definitions  one  finds 


that  oT 


p  2  ,  2 

cTp  =  1/LC  -  a  and  osL 


(i  +  p2)"1/2(L/C)1/2;  thus 


m 


/  i  /t  \l/^ 

v •-L; 


r  -1+  -1 

-  p  tan  p  - 


(10) 


The  square  of  this  expression  corresponds  to  that  given  by  Anderson  and  Smith 

except  that  their  exponential  term  is  incorrect  and  an  incorrect  factor  of 

t/5  enters  because  of  their  failure  to  evaluate  at  t  .  Maximum  power  is 
'  m 


obtained  when  m  =  0,  f  =  1.10  and  p 
P 


Pi 


Thus 


max 


1.10  V2  (C/L r'2  e-1’50 


(U) 


For  a  typical  case,  L  =  l/3  nh.,  C  =  0.3  iif,,  V  *  28  k  v  and  as  before 

/  \l/2 

Pt  ■  1.515*  With  these  values  (L/C)  '  -  1.05  ohms,  i  a  15.3  k  a  and 

P  *  2.24  x  10°  watts.  If  approximately  this  power  is  consumed  during 
l3*  sec.,  around  the  peak,  then  22  Joules  are  used  during  the  interval. 

With  the  exploding  wire  replaced  by  a  copper  strip,  the  circuit  has  a 

logarithmic  decrement  given  approximately  by  5  *  ln2  *  2n/p  from  which 

computation  gives  R  =  0.23  ohms.  For  maximum  power  at  current  peak, 

T_  /o 

R  =  1.10  (L/C)J'/ *"  »  1.1$.  .  Thus  for  optimum  conversion  of  the  condenser 

energy  in  the  explosion,  the  wire  should  furnish  about  0.9  ohms. 
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2 ;  INTEGRATED  POWER. 

Since  available  evidence  indicates  that  the  circuit  containing  an 
exploding  wire  of  nearly  optimum  size  is  broken  sometime  during  the  latter 
part  of  the  first  half  cycle,  it  is  of  interest  to  inquire  whether  power 
consumed  during  that  interval  can  be  maximized  by  selection  of  an  appropriate 
circuit  resistance.  To  this  end  the  energy  taken  from  the  condenser  at  time 


W  = 


.  m 
/o' 


i2  R  dt 


This  equation  can  be  integrated  in  compact  form  by  letting  f/2 

=1  to 

and  m  =  (10)  JV‘_j  thus 


(12) 
g  =  sin0 


R  ■  2  sin  0  (L/C)1/2 


n  4  m  ^ 

bxu  y 


and 


0)  =  cu  cos  0 
o 

p  “  to/a  -  cot  0 


(12a) 

/  i  nv  \ 

\XC10) 

(12c) 

(I2d) 


U ^  +\Vl  t.VlOfSO  A 1  Q  ona  Ptl/1  a nmo  ol  orriftn  +  owr  cfr^n/^n1fi+.v^  nol  +  s-M-ta 

"  *  “*•  w  A  WA,V«V4kV**W  M  WM1V  VO,  VU1V44  VMA  j  WA  AQV*4V«aV  WA  AWWA  Wi  UUOl  VI  iUU  VI  WUO 


-n 


(1  ’d)  integrates  to 

W  «  W q  ^1  -  sec^0  e“cpXT'anJw  +  si n0  sin(2a>T  -  0)^  ,  (lj) 

where  WQ  (■  CV2/2)  represents  the  initial  energy  in  the  condenser.  From 
the  earlier  requirement  on  f  note  that  0  sin0  ^ lj  thus  the  negative 

O 

term  in  (13)  always  remains  negative.  From  the  fact  that  the  power  i'lR  is 
positive  definite  it  follows  that  W  is  a  non-decreasing  function  of  T. 

To  examine  W  for  stationary  points  with  respect  to  T  and  R,  employ 

ft  =«ln  |(W  -  W)/W  "1  for  convenience.  With  this  definition  the  first 

*-  °  °-*  a  t  /o  a 

tvc^oId  hkb  PttJ/tVP  _  U  ( ~>.C\  \  iTJ  /in  -  .  W  l  ),T  In  \ ^  J* 

*-%*.  w  f  \J  J.  —  —  n  V-  \  UH/  W-*-  /  cajliu*.  U»'/  ”  n  \  *+J_y  y  OCVJJfc/  C  \  Q)£/  }  • 

ft  ^  °  °  Q 

since  e'/yQ,  the  stationary  points  will  be  given  by  setting  e"*  dft/dT 
»  see  0  e^  OQ/d0)  =  0.  From  the  fact  that  e*^— >0  as  T-**>it  is  clear  that 
as  time  increases  the  W  surface  approaches  the  plane  W  =  W  for  any  value 


Performing  the  partial  differentiation  leads  to  two  equations 


sin  0  (cos  2  coT  -  1)  =  0 


and 


sec  0  j  sin  2(cuT-0)  +  2tan0  |l  +  sin  0  sin.’  (£h)T-0)J 

A  m  ft  j -  rl  >  A  itTI*]  a 

-  aax  |^l  +  uanp  sin  £iDxj^  =  u 
A  number  of  cases  arise  from  these. 


(14) 


(15; 
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In  the  first  place  (ll+)  indicates  that  stationary  points  in  time  occur 

When  dl  ■  rut.  that  is  at  points  where  the  current  is  zero.  It  follows  that 

the  curves  of  W  vs  T  have  horizontal  tangents  at  these  points,  beginning 

with  T  *  0  for  n  =  Q,  and  rise  smoothly  in  a  series  of  steps  to  the  maximum 

W  .  The  size  of  the  steps  is  not  constant  but  diminishes  montonically  as 
o 

integer  n  increases. 


When  sin  0  =  0  (R  *  0),  Eq.  (14)  is  satisfied  and  (15)  leads  to  the 


r\Y\  onoo  *"P 

\y  wii  .  " 


V  U14U  J. /  —  V  “  \J  f  y  UUC 


origin  is  clearly  a  true  minimum. 


T  -f* 


,4*  (k  1  C\  +  Viori  oa  /ivHP 

7X11  p  J~  \y  y  WllV-il  uo  wux  vx  ^  U7X 


Imi+  -I, 


r\r\  "1  ir  C’O-f-H  a  -fM  a  A  ^  -f* 

V11XJ  OUUXEUXVU  JLX 


n  ■  0.  Thus  the  equation  off  ■  aiQ  cos  jz>  T  =  0  must  hold.  With  cos  jz>  f  0, 
the  0  axis  emerges  as  a  minimum  line  of  the  surface.  This  is  to  be  expected 
«1nne  no  nower  has  been  consumed  at  T  =  0  for  anv  value  of  0  ( resistance^  . 


mi _ _ li _  _i _ * _ _ a _ j _ _ _ d(  m  _  a  j  _  x _ _ x.  ^ 

me  uuxier  exiuj.ee  rexiu.ejrj.xig  eue  yj  jl  =  u  jlb  uu  ecu  // 


—  /  O  J  _  _1 _ _ _ XI _ 

life.  i»c,  euuuee  uxie 


critical  damping  resistance.  But  this  value  fails  to  satisfy  (15)  because 
sec  0-^ooas  0-*  it/2.  Investigation  of  the  limit  shows  that  the  critical 


rm-i  -4  ri  rr  ahvivo  ^  i 
lupxiip  wmx  » v  xi 


no'b  s  stationary  iins  of  tlis  W  surfscs , 


Our  analysis  shows  that  W  has  no  true  maxima  or  minima  for  Interior 


iraliiaa  r\-P  +Vin  raa<1  otonnfl  • 
vaxu^o  ui  uuv.  x  ox  a  voiivv  y 

the  energy  consumed  by  some  finite  tiie.  On  the  other  hand  for  some  preferred 
value  of  T,  the  W  curve  may  have  a  relative  maximum  as  a  function  of  0. 

To  investigate  this  possibility,  we  look  for  solutions  to  ( 15 )  other  than 
the  one  already  found. 

.QiiTmrtaa  -f'r'.v*  ovcntml.  ito  cold/tf.  fhftcp  fnr  r»Vi  /I'/P  —  ( Pm  a.  1  tr  /p 

kj  u^vuv.  y  x  vx  vamiu^x.v  y  nw  vvx.wv  w  »-«  *  w*  ”  v»**  w-*-  —  \  «  ■  ■*■/“/  •— 


0,  1, 


2...  o  If  in  addition  0  <  n/2  is  required,  then  the  values 


with  m 

of  T  so  defined  are  finite  and  correspond,  according  to  (l),  to  the  points 
of  tanflrenev  of  the  current  curve  with  its  exnonentiailv  decreasing  envelopes 


xnese  points  occur  later  in  uime  uxitui  une  cun’etsponuiug  curreuu  mtunma,  or 
minima,  i.e„,  later  in  time  them  the  power  maxima,  and  thus  represent  em 
integration  to  times  well  past  the  points  of  maximum  power.  The  shortest 


4  - -  1  m  - V.,, 

ilibci  vai  X  y  ftXVCU  UJ 


_  r\  *,4  11  4-  4  „4- ~4- 

—  \j  y  will  uc  uxic  iuwou  inucicooiug  iium  o  ucuiupuixi  u 


of  the  exploding  wire.  It  is  worth  noting  that  for  ft  =  Jt/2,  the  angular 
frequency  vanishes,  i.e.,  co  =  0  and  all  T  values  for  which  dT  =  const,  as 
0  =xjt/2  coincide  with  the  point  at  infinity. 
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For  the  points  of  tangency,  (15)  vanishes  if  tan  0  *  (2  m  +  l)«/U.  When 
m  =  0  we  find  sin  0  =  .618  and  R  =  l^UCL/C)1^.  Since  this  value  is  about 
10$  larger  than  the  value  of  R  which  gives  maximum  power  consumption  at 
peak  current,  we  suspect  a  relative  maximum  of  energy  consumption,  i.e., 
a  maximum  in  the  W  vs  R  curve.  Examination  of  (13)  for  values  in  the 
neighborhood  of  0  =  tan-1(n/4)  while  cdT  =  n/2  shows  that  this  is  indeed  the 
case.  Thus  if  integrated  power,  in  other  words  total  energy,  released  in 
the  wire  is  the  criterion  of  importance;  then  the  optimum  resistance  R  for 
the  first  point  of  tangency  is  slightly  more  than  10$  larger  than  the  value 
producing  maximum  power  release  at  peak  current.  If  the  wire  should  last 
longer  than  the  first  half  cycle,  then  even  larger  values  of  R  would  be 
appropriate  corresponding  to  m  *  1;  and  in  the  limit  of  indefinitely  many 


rvvrr»"l  o  c 
V  j  L 


15  or\r»w\OrtVioci  4*  A  /■»  rs  1  rJowm-1  wrr  irnl  hq 

xv  vuv.  jl.  umuy  v  >.  u.v^  • 


Should  other  values  of  Interval  T  be  chosen,  then  it  is  likely  that 
different  optimum  values  for  R  would  be  found,  and  it  seems  probable  that 
these  values  will  be  larger  or  smaller  than  those  Just  found  depending  on 
whether  o>T  is  taken  larger  or  smellier  than  (2  m  +  l)*/2. 


3.  DISCUSSION 


The  conclusions  which  emerge  from  the  preceding  analysis  of  the  under - 
damped  R-L-C  circuit  are  as  follows: 

(1)  For  maximal  power  dissipation  at  .the  first  current  peak  a  circuit 
resistance  R  »  l.lOfL/C)1^  should  be  used.  This  more  precise  value  is  10$ 
higher  than  that  given  by  Anderson  and  Smith  in  1926. 

(2)  For  maximum  energy  dissipation  in  the  interval  terminated  by  the 

first  contact  of  the  current  curve  and  its  exponential  envelope,  the  re- 

l/2 

oistance  should  be  R  =  l.2k(L/C)  '  .  As  the  interval  increases,  the  optimum 
value  of  resistance  will  increase  toward  the  critical  damping  resistance. 

It  has  not  been  our  purpose  here  to  demonstrate  the  relevance  of  these 
criteria  for  the  experimental,  exploding  wire  phenomenon.  The  complexity 
of  this  process  flefles  exact  analysis  at  present.  This  much  may  be  said; 

i 

the  wire  resistance  is  certainly  not  constant  during  the  explosion,,,  but 
must  change  rapidly  and,  in  all  probability,  discontinuous ly  with  time. 

It  may  be  conjectured  that  if  the  process  is  rapid  enough  to  be  essentially 
adi  aibaviC  uuSIl  the  larger  value  of  (2)  will  apply;  while  if  the  explosion 
depends  mainly  on  heat  generated  near  the  first  current  peak  then  the  smaller 
value  of  (l)  is  appropriate.  In  intermediate  cases  the  inequality  1.10 
±  rCl/G)"1' 2  ■£;  1.2k  may  be  satisfied. 
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Ballistic  Research  laboratories,  APG 

ENERGY  PARTITION  IN  THE  EXPLODING.  WIRE  PHENOMENON 

F.  JO.  Bennett 

BRL  Report  No.  1056  October  19581 

DA  ;?roj  No.  5B03--03-001,  ORD  ProJ  TB3-O1O0 
UNCLASSIFIED  Report 

Streak  cameraaid  oscilloftraphic  circuit  dancing  data  Sire  presented  for 
copjjer  wires  varying  in  diameter  from  3  to  8  mils.  A  maximum  of  specific:  shock 
wave  energy  in  the  Induced  flow  is  found  at  a  wire  diameter  different  from  that 
of  a  minimum  in  the  total  damping;  time  of  the  circuit.  'This  displacement  is 
shown  to  be  caused  by  the  presence  of  residual  circuit  resistance.  Tide  argument 
is  based  on  a  critical  analysis  of  optimum  damping  conditions  in  the  (exploding 
wire  circuit.  A  maximum  of  apparent  energy  within  the  contact  surface  appears 
at  about  the  same  wire  diameter  as  the  minimum  of  total  dangling  time.  Discussion 
of  the  implications  of  the  Taylor- Lin  similarity  theory  indicates  that  lack  of 
similarity  of  the  flow  is  probably  connected  with  the  displacement  of  the  maximum 
energies  associated  with  shock  ws.ve  and  contact  surface. 
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DA  ProJ  No.  5B03-03-001,  ORD  ProJ  TB3-0108 
UNCLASSIFIED  Report 

Streak  camera  end  oscillographic  circuit  lamping  data  are  presented  for 
copper  wires  varying  in  diameter  frem  3  to  8  mils.  A  maximum  of  specific  shock 
wave  energy  in  the  induced  flew  is  found  at  a  wire  diameter  different  from  that 
of  a  minimum  in  the  total  damping;  time  of  the  circuit.  This  displacement  is 
shown  to  be  caused  by  the  presence  of  residual  circuit  resistance.  The  argument 
is  based  on  a  critical  analysis  of  optimum  damping  conditions  in  the  exploding 
wire  circuit.  A  maximum  of  apparent  energy  within  the  contact  surface  appears 
at  about  the  same  wire  diameter  as  the  minimum  of  total  damping  time.  Discussion 
of  the  implications  of  the  Taylor-Lln  similarity  theory  indicates  that  lack  of 
similarity  of  the  flow  is  probably  connected  with  the  displacement  of  the 
maximum  energies  associated  with  shock  wave  and  contact  surface. 
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UNCLASSIFIED  Report 

Streak  camera  and  oscillograptilc  circuit  dashing  data  are  presented  for 
copper  wires  varying  In  diameter  from  3  to  8  mils.  A  maximum  of  specific 
shock  waire  energy  in  the  Induced  flow  is  found  at  a  wire  diameter  different 
from  that  of  a  minimum  in,  the  total,  damping  time  of  the  circuit.  This  displace¬ 
ment  is  shown  to  be  caused  by  the  presence  of  residual  circuit  resistance. 

The  argument,  is  based  on  a  critical,  analysis  of  optimum  damping  conditions 
In  the  exploding  wire  circuit.  A  maximum  of  apparent  energy  within  the  contact 
surface  appears  at  about  the  same  wire  diameter  as  the  minimum  of  total  damping 
time.  Discussion  of  the  iipli  cations  of  the  Taylor-Lln  similarity  theory 
indicates  ttiat  lack  of  similarity  cf  the  flow  is  probably  connected  with  the 
displacement  of  the  maxiann  energies  associated  with  shock  wave  and  contact 
surface. 
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Streak  camera  ancl  osciliograpMc  circuit  damping  data  are  presented  far 
copper  wires  varying  In  diameter  from  3  to  8  mils.  A  maximum  of  specific 
shock  wave  energy  in  the  induced  fl.ow  is  found  at  a  wire  diameter  different 
from  that  of  a  minimum  in.  the  total,  dating  time  of  the  circuit.  Tills  dis¬ 
placement  is  shown  to  he  caused  by  the  presence  of  residual  circuit  resistance. 
The  argument  is  based  on  a  critical,  analysis  of  optimum  damping  conditions 
in  the  exploding  wire  circuit.  A  naxlaus  of  apparent  energy  within  the  contact 
surface  appears  at  about  the  same  wire  diameter  as  the  minimum  of  total  damping 
time.  Discussion  of  the  implications  of  the  Taylor-Lln  similarity  theory 
indicates  thiat  lack  of  similarity  cf  the  flow  is  probably  connected  with  the 
displacement,  cf  the  maximum  'energies  associated  with  shock  wave  and  contact 
surface. 


